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Abstract 

A general class of loop quantizations for anisotropic models is introduced and 
discussed, which enhances loop quantum cosmology by relevant features seen in in- 
homogeneous situations. The main new effect is an underlying lattice which is being 
refined during dynamical changes of the volume. In general, this leads to a new fea- 
ture of dynamical difference equations which may not have constant step-size, posing 
new mathematical problems. It is discussed how such models can be evaluated and 
what lattice refinements imply for semiclassical behavior. Two detailed examples 
illustrate that stability conditions can put strong constraints on suitable refinement 
models, even in the absence of a fundamental Hamiltonian which defines changes of 
the underlying lattice. Thus, a large class of consistency tests of loop quantum grav- 
ity becomes available. In this context, it will also be seen that quantum corrections 
due to inverse powers of metric components in a constraint are much larger than they 
appeared recently in more special treatments of isotropic, free scalar models where 
they were artificially suppressed. 



1 Introduction 

Loop quantum cosmology [Ij was designed to test characteristic effects expected in the full 
framework of loop quantum gravity Pi El H] . Implementing symmetries at the kinematical 
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quantum level allows explicit treatments of the dynamical equations while preserving basic 
features such as the discreteness of spatial geometry [5]. (See also [SI [TJ [H [H [TOl [11] for 
recent work on symmetry reduction in quantum theories.) Indeed, several new, initially 
surprising results were derived in different applications in cosmology and black hole physics. 
By now many such models have been studied in detail. 

As the relation of dynamics to that of a possible full framework without symmetries is 
not fully worked out, detailed studies can be used to suggest improvements of the equa- 
tions for physically viable behavior. Comparing results with full candidates for quantum 
dynamics can then provide stringent self-consistency tests of the overall framework. It is to 
be seen if, and how, such alterations of quantization procedures naturally result from a full 
quantization. The first example of this type related to the stability behavior of solutions to 
the difference equations of isotropic loop quantum cosmology, which was studied in [121 [13] 
and was already restrictive for models with non-zero intrinsic curvature. Another limita- 
tion, realized early on [2], occurs in the presence of a positive cosmological constant A. 
In an exact isotropic model, the extrinsic curvature scale is given by /c = d = A/STrGa^A/S 
which, due to the factor of a^, can be large in a late universe although the local curva- 
ture scale A might be small. Extrinsic curvature plays an important role since in a flat 
isotropic model it appears in holonomies on which loop quantizations are based in such a 
way that only e*"^ with a G M can be represented as operators, but not k itself [15] . Large 
values of k would either require one to use extremely small a in the relevant operators, or 
imply unexpected deviations from classical behavior. In fact, holonomies as basic objects 
imply that the Hamiltonian constraint is quantized to a difference rather than differential 
equation p!6j since k in the Hamiltonian constraint (as in the Friedmann equation) is not 
directly quantized but only exponentials e*°'^. These are shift operators instead of differen- 
tial operators. For a large, semiclassical universe a Wheeler-DeWitt wave function should 
be a good approximation to the basic difference equation of loop quantum cosmology [17] 
which, in a representation as a function of the momentum p = conjugate to fc, would be 
oscillating on scales of the order (oa/A)"^. This scale becomes shorter and shorter in an 
expanding universe, eventually falling below the discreteness scale of the difference equa- 
tion of loop quantum cosmology. At such a point, discreteness of spatial geometry would 
become noticeable in the behavior of the wave function (independently of how physical 
observables are computed from it) although the universe should be classical. 

This does not pose a problem for the general formalism, because it only shows that the 
specific quantization of the exact isotropic model used reaches its limits. Physically, this 
can be understood as a consequence of a fixed spatial lattice being used throughout the 
whole universe evolution. Exponentials e*"'^ in isotropic models derive from holonomies 
hf.{A) = V exp{J^Al^Tie°'dt) of the Ashtekar connection along spatial curves e. All the 
freedom contained in choosing edges to capture independent degrees of freedom of the full 
theory reduces, in isotropic models, to the single parameter a which suffices to separate 
isotropic connections through all functions e*"'^. The parameter a, from the full perspective, 
is thus related to the edge length used in holonomies. Using a fixed and constant a is 
analogous to using only edges of a given coordinate length, as they occur, for instance, in 
a regular lattice. In the presence of a positive cosmological constant, for any a a value of k 
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will then be reached such that e differs strongly from iak. From the lattice perspective, 
this means that the local curvature radius becomes comparable to or smaller than the 
fixed lattice scale corresponding to a. Such a fixed lattice ceases to be able to support all 
small-scale oscillations relevant for a semiclassical geometry. 

This is not problematic if it occurs in a quantum regime where dynamics is indeed 
expected to differ from the classical one, but it poses a problem in semiclassical regimes. 
A better treatment has to refer to changing lattices, which is not easy to implement in 
a straightforward quantization of purely homogeneous models. In a dynamical equation 
closer to what is expected from the full framework, lattice refinements would take place 
during the evolution since full Hamiltonian constraint operators generally create new ver- 
tices of a lattice state in addition to changing their edge labels [HI [19]. While k increases 
with increasing volume, the corresponding a decreases since the lattice is being refined all 
the time. For a suitable lattice refinement, the increase in k can be balanced by the decrease 
of a such that ak stays small and semiclassical behavior is realized for any macroscopic 
volume even with A > 0. This provides an interesting relation between the fundamental 
Hamiltonian, which is responsible for the lattice refinement, and semiclassical properties 
of models. Testing whether an appropriate balance between increasing k and lattice re- 
finements can be reached generically can thus provide stringent tests on the fundamental 
dynamics even without using a precise full Hamiltonian constraint operator. 

This feature of lattice refinements was not mimicked in the first formulations of loop 
quantum cosmology [Sni EH El ESI [IS] since the main focus was to understand small- volume 
effects such as classical singularities [231 (21] • In this context, lattice refinements appear 
irrelevant because only a few action steps of the Hamiltonian, rather than long evolution, 
are sufficient to probe a singularity. By now, perturbative regimes around isotropic models 
have been formulated in loop quantum cosmology which are inhomogeneous and thus must 
take into account lattice states and, at least at an effective level, lattice refinements [H]- 
One special version, corresponding to lattices with a number of vertices growing linearly 
with volume in a specific way referring to the area operator, has been studied in detail in 
isotropic models with a free, massless scalar [25]. Although the complicated relation to a 
full, graph-changing Hamiltonian constraint is still not fully formulated, such models allow 
crucial tests of the local dynamics. 

While isotropic models can easily be understood in terms of wave functions on a 1- 
dimensional discrete minisuperspace in terms of oscillation lengths [26] , anisotropic models 
with higher-dimensional minisuperspaces can be more subtle. In such models, limitations 
similar to that of a cosmological constant have been observed as possible instabilities 
of solutions in classical regions or the lack of a sufficient number of semiclassical states 
[271 [2H1 [29]. For the partial difference equations of anisotropic models in loop quantum 
cosmology, stability issues can be much more severe than in isotropic models and thus 
lead to further consistency tests which might help to restrict possible quantization freedom 
(see, e.g., [30]). In this paper we therefore introduce the general setting of anisotropic 
models taking into account lattice refinements of Hamiltonian constraint operators, fo- 
cusing mainly on the anisotropic model which corresponds to the Schwarzschild interior. 
As we will see, the type of difference equations in general changes since they can become 
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non-equidistant. This leads to new mathematical problems which we address here briefly, 
leaving further analysis for future work. The examples presented here already show that 
one can distinguish different refinement models by their stability properties. The refine- 
ment model corresponding to [25j turns out to give unstable evolution of the Schwarzschild 
interior, while a new version, whose vertex number also grows linearly with volume, is 
stable. Compared to isotropic models which are sensitive only to how the vertex number 
of a state changes with volume, anisotropic models allow one to test much more detailed 
properties. 

An appendix discusses subtleties in how homogeneous models faithfully represent inho- 
mogeneous states, mainly regarding the magnitude of corrections arising from quantizations 
of inverse metric components which often plays a large role in cosmological applications. 

2 Difference equation for the Schwarzschild interior 
with varying discreteness scale 

Basic variables of a loop quantization are holonomies along lattice links and fluxes over 
transversal surfaces. For the Schwarzschild interior [31], the connection used for holonomies 
and the densitized triad used for fluxes take the form 

A^Tjdx" = CT^dx + (ari + bT2)d{} + {—bri + 0x2) sinMcy^ + cos-i^dc/? (1) 
■ d _ i9__ d __i9 

E^t'—- = PcT3Smi^— + {paTi + PbT2) sin^9— + {-pbTl+PaT2)Tr ■ (2) 

OX'^ OX ov oip 

Coordinates (x, i), ip) are adapted to the symmetry, with polar angles and cp along orbits 
of the rotational symmetry subgroup, and tj = —^crj in terms of Pauli matrices. Spatial 
geometry is determined by the spatial line element, which in terms of the densitized triad 
components is 

ds2 = 4^dx2+|p,|dfi2 (3) 

\Pc\ 

obtained from g*^* = E^E\/\ deti?||. We will also use the co-triad e^, i.e. the inverse of 
et = E-/^\detE% 

elTidx" = ecTgdx + {can + ebT2)d{} + {-CbTi + CaTa) sinMip (4) 
with components 

_ sgnpc VPl + Pb _ \lW\Ph ^ ^ -JW^Pa 



\fWl ' \/p1^pI \/pI + Ph 

The phase space is spanned by the spatial constants {d,b,c,pa,Pb,Pc) G IR^ with non- 
vanishing Poisson brackets 

{d,pa} = 'jG/Lo , {b, pb} = -fG / Lq , {c,pc} = 2'jG/Lo 
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where G is the gravitational constant and 7 the Barbero-Immirzi parameter 
Moreover, Lq is the size of a coordinate box along x used in integrating out the fields in 



87r7G J " ' 276- " 7G " 76- 

to derive the symplectic structure. The SU(2)-gauge transformations rotating a general 
triad are partially fixed to U(l) by demanding the x-component of to point in the 
internal rs-direction in ([2]). The U(l)-gauge freedom allows one to set a = = p^, still 
leaving a discrete residual gauge freedom {b,ph) ^ (—6, —pb)- The remaining variables can 
be rescaled by 

{b, c) := {b, Loc) , {pb,Pc) ■-= (LoPb^Pc) ■ (6) 
to make the canonical structure Lg-independent: 

{b,pb}=iG , {c,pJ = 27G. (7) 

This rescaling is suggested naturally by holonomies, as written below, and fluxes which are 
considered the basic objects in loop quantizations. 

To express the elementary variables through holonomies, which unlike connection com- 
ponents will be promoted to operators, it suffices to choose curves along the x-direction of 
coordinate length tLq and along ^ of coordinate length /x since this captures all information 
in the two connection components, 

(A) = exp / dxcTs = cos h 2t^ sin — (8) 

Jo 2 2 

hi^^ (A) = exp r di&bT2 = cos ^ + 2t2 sin ^ . (9) 
Jo 2 2 

The quantum Hilbert space is then based on cylindrical states depending on the connection 
through countably many holonomies, which can always be written as almost periodic func- 
tions c) = T ffj;-r ^^P f (/^^ + '^^) variables. These form the set of functions 
on the double product of the Bohr compactification of the real line, which is a compact 
Abelian group. Its Haar measure defines the inner product of the (non-separable) Hilbert 
space, in which states 

(6, c|/i, r) = (^''+"") fi,TeR. (10) 

form an orthonormal basis. Holonomies simply act by multiplication on these states, while 
densitized triad components become derivative operators 

^^ = -^74^ Pe = -2^74^ (11) 

using the Planck length £p = \J Gh. They act as 

Pfe|/i,r) = i74/i|/i,r), = 74T|/i,r) , (12) 
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immediately showing their eigenvalues. 

To formulate the dynamical equation, one has to quantize the Hamiltonian constraint 

1 r rpai rpbj 

H=- d'xe,,,{-Ft, + ^'Q>:,)-=== (13) 
7 J ^|detE| 

where fi^j^r^da;" A dx'' = — sint^rsdt? A dcp is the intrinsic curvature of 2-spheres, while 
is the curvature computed from ignoring the spin connection term sin-^r^dip. Following 
standard procedures a Hamiltonian constraint operator can be expressed in the basic op- 
erators. First, one replaces the inverse determinant of by a Poisson bracket, following 

= E i^""-f4^{4^-\n (14) 



^|detE| 4:TT^G 



K&{x,'d,ip} ° 



with edge lengths £g = SLq and Iq^'^ = 6, and left-invariant 1-forms on the symmetry 
group manifold. For curvature components one uses a holonomy around a closed loop 



Fl,(x)r, = ^{hfj - 1) + 0((6^ + cy/^VA) (15) 

with 

i,fj = h'Phf(h'P)-\hfr' (16) 

and Aij being the coordinate area of the loop, using the corresponding combinations of 
£q. In these expressions, a parameter 5 has been chosen which specifies the length of edges 
with respect to the background geometry provided by the symmetry group. Putting all 
factors together and replacing Poisson brackets by commutators, one has 

\IJK 

A-/ 3^3n2\-lfa ■ Sb . 6c Sc f . Sb^y 6b 6b . 6b 

= 4z(7 I 8 sm — cos — sm — cos — sm — v cos cos — V sm — 

^'P^V 2 22 2V2 2 2 2 

f,.2^b o6b 2c2\ ( ■ (^c- 6c 6c^y . 6cW 

+ I 4 sm — cos h 7 I I sm — V cos cos — V sm — II (17) 

V 2 2 yy2 2 2 ^ J J 

which acts as 

H^'^\I2,t) = {2^'6Hl)-'[2{V^+s,r-V^-s,r) (18) 
X (l/i + 2(5, r + 26) -\fi + 26,T- 26) - - 25, r + 26) + |/i - 25, r - 26)) 
+ {V^,r+s - V^,r-5){\^i + 45, r) - 2(1 + 27^52) r) + |/x - 45, r))] 

on basis states. This operator can be ordered symmetrically, defining ^iymm := \{H''^'' + 
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_f)'('^)t)^ whose action 

-(V^+5,r - V'^-5,T + V'm+35,t-25 - V^+s,T-2s)\f^ + 2(5, T - 2(5) 
-{V^+5,T - V^-S,T + Vf,^s,T+25 - V'^-35,T+25)|/i - 2(5, T + 2(5) 
+ (V;.+5,T - V'/,_5,T + V'm-5,t-25 - Vf,^3S,T~2s)\fJ' - 2(5, T - 2(5) 
+ |(V^,T+5 - ViJ.,T-S + V'^+45,r+5 - V'^+45,T-5)|/i + 45, r) 

-2(1 + 27^(52) (\/^,.+5 - r) 

+ |(V^,r+5 - V^^r-S + V;_45,r+5 " Vf,.4S,r-s)\fi " 4(5, t)] . (19) 

Transforming this operator to the triad representation obtained as coefficients of a wave 
function = J^/xt '^tJ.Al^i^) ™- the triad eigenbasis and using the volume eigenvalues 



V^^r = 47ry/|(pJ^,.|(p,)^,. = 27r(74)3/2^;,, 

a difference equation 

^3/2^3 

^^(4?LlV'))M,r = 2(5(v/|r + 2(5| + V^) {'iP^.+25,r+2S - V^M-25,r+25) 

+ W\r + 5\ - ^\t-5\) ((/i + 25)^^+45,r - 2(1 + 272(52)/i^^,, + (;x - 2(5)^^_45,0 
+2(5(v^7^+ V^) (^^ 
= (20) 

results for physical states. (For small the equation has to be specialized further due to 
the remaining gauge freedom; see |3]]|. This is not relevant for our purposes.) 



2.1 Relation to fixed lattices 

Although there are no spatial lattices appearing in the exactly homogeneous context fol- 
lowed here, the construction of the Hamiltonian constraint mimics that of the full theory. 
States are then associated with spatial lattices, and holonomies refer to embedded edges 
and loops. The parameter 5 is the remnant of the loop size (in coordinates) used to act 
with holonomies on a spatial lattice. As one can see, this parameter is important for the 
resulting difference equation, determining its step-size. The above construction, using a 
constant 5, can be seen as corresponding to a lattice chosen once and for all such that the 
loop size is not being adjusted even while the total volume increases. As described in the 

^Note that the first factor of 2 in the next-to-last fine was missing in and analogous places in 
subsequent formulas. This turns out to be crucial for the stability analysis below. In particular, with the 
corrected coefficient the quantization of the Schwarzschild interior in [31] is unstable for all values of 7. 
Possible restrictions on 7, as suggested in ^3(T based on a difference equation with the wrong coefficient, 
then do not follow easily but could be obtained from a more detailed analysis. 
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introduction, this ignores the possible creation of new lattice vertices and links, and can 
be too rigid in certain semiclassical regimes. 

To express this clearly, we now construct holonomies which are not simply along a single 
edge of a certain length 5, but which are understood as holonomies along lattice links. We 
keep our coordinate box of size Lq in the x-direction as well as the edge length io. If this 
is a link in a uniform lattice, there are j\4 — Lq/Iq lattice links in this direction, and a link 
holonomy appears in the form 

= exp(£ocr3) = cxp^ioCTs/ Lq) = cxp{cTs/M^) (21) 

when computed along whole lattice edges. Thus, a constant coefficient in holonomies 

corresponds to a fixed lattice whose number of vertices does not change when the volume 
increases. Lattice refinements of an inhomogeneous lattice state, on the other hand, can be 
mimicked by a parameter which depends on the phase space variables, most importantly 
the triad components. If this is carried through, as we will see explicitly below, the step-size 
of the resulting difference equation is not constant in the triad variables anymore. 

2.2 Lattice refinements 

Let us now assume that we have a lattice with A/" vertices in a form adapted to the 
symmetry, i.e. there are vertices along the x-direction (whose triad component pc gives 
rise to the label r) and J\f^ vertices in spherical orbits of the symmetry group (whose triad 
component pt, gives rise to the label //). Thus, jV = NxM^. 

Since holonomies in such a lattice setting are computed along single links, rather than 
through all of space (or the whole cell of size Lq), basic ones are = cxp(£Qcr3) and 
h{t = exp(£Q6r2), denoting the edge lengths by and keeping them independent of each 
other in this anisotropic setting. Edge lengths are related to the number of vertices in 
each direction by £g = Lq/Mx and ^Q = I/M-q. With the rescaled connection components 
c — Lqc and b — b we have basic holonomies 

hx = exp(£oLo ^cTs) = exp(cT3/A4) , = exp(£o6T2) = exp(6T2/M») . (22) 

Using this in the Hamiltonian constraint operator then gives a difference equation whose 
step-sizes are l/A/}. 

So far, we only reinterpreted 6 in terms of vertex numbers. We now turn our attention to 
solutions to the Hamiltonian constraint which, in the full theory, usually changes the lattice 
by adding new edges and vertices while triad eigenvalues increase. For larger /i and r, the 
Hamiltonian constraint thus acts on a finer lattice than for small values, and the parameter 
A/" for holonomies appearing in the constraint operator is not constant on phase space but 
triad dependent. Due to the irregular nature of lattices with newly created vertices such a 
refinement procedure is difficult to construct explicitly. But it is already insightful to use 
an effective implementation, using the derivation of the Hamiltonian constraint for a fixed 
lattice, but assuming the vertex number A/'(/i, r) to be phase space dependent. Moreover, 
we include a parameter S as before, which now takes a value < S < 1 and arises because 
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a graph changing Hamiltonian does not use whole lattice edges but only a fraction, given 
by 5B Effectively assuming in this way that the lattice size is growing through the basic 
action of the Hamiltonian constraint, we will obtain a difference equation whose step-size 
6/// is not constant in the original triad variables. 

For the Schwarzschild interior, we have step sizes S/N"^ for fi and 6/J\fx for r. Go- 
ing through the same procedure as before, we end up with an operator containing flux- 
dependent holonomies instead of basic ones, e.g., A/l.(/i, r)/ij. = A/l.(/i, r) exp(cr3/A4(/i, r)) 
which reduces to an A4-independent connection component c in regimes where curvature 
is small. Keeping track of all prefactors and holonomies in the commutator as well as the 
closed loop, one obtains the difference equation 



(A total factor MxM^ for the number of vertices drops out because the right hand side is 
zero in vacuum, but would multiply the left hand side in the presence of a matter term.) 

3 Specific refinement models 

For further analysis one has to make additional assumptions on how exactly the lattice 
spacing is changing with changing scales and r. To fix this in general, one would have 
to use a full Hamiltonian constraint and determine how its action balances the creation 
of new vertices with increasing volume. Instead of doing this, we will focus here on two 
geometrically motivated cases. Technically simplest is a quantization where the number of 
vertices in a given direction is proportional to the geometrical area of a transversal surface. 
Moreover, the appearance of transversal surface areas is suggested by the action of the full 
Hamiltonian constraint which, when acting with an edge holonomy, creates a new vertex 
along this edge (changing A/} for this direction) and changes the spin of the edge (changing 

precise value can be determined only if a precise implementation of the symmetry for a fixed full 
constraint operator is developed. Currently, both the symmetry reduction for composite operators and a 
unique full constraint operator are lacking to complete this program and we have to work with 5 as a free 
parameter. This parameter is sometimes related to the lowest non-zero eigenvalue of the full area operator 
jTSl I25j. From the inhomogeneous perspective of lattice states used here, however, there is no indication 
for such a relation. 




+ (/i - 26J\f^{fi, t) ^)^^-45Ar^(^,r)-l,r) 

0. 



(23) 



with 




(24) 
(25) 
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the area of a transversal surface). It also agrees with [251 [3lj, although the motivation in 
those papers, proposing to use geometrical areas rather than coordinate areas Au in (flGll . 
is different. 

Geometrically more intuitive is the case where the number of vertices in a given direction 
is proportional to the geometrical extension of this directionjfl The resulting difference 
equation will be more difficult to deal with due to its non-constant step-size, but naturally 
gives rise to Misner-type variables. This case will also be seen to have improved stability 
properties compared to the first one using areas. In both cases, M cc V is assumed, i.e. 
the lattice size increases proportionally to volume. This is not necessary in general, and 
we choose these two cases mainly for illustrative purposes. In fact, constant Af as in [15] 
and N' oc V first used in [25] are two limiting cases from the full point of view, the first 
one without creating new vertices and the second one without changing spin labels along 
edges since local lattice volumes V/A/ remain constant. In general, both spin changes and 
the creation of new vertices happen when acting with a Hamiltonian constraint operator. 
Thus, one expects J\f oc V°' with some < a < 1 to be determined by a detailed analysis 
of the full constraint and its reduction to a homogeneous model. Even assuming a certain 
behavior of J^iV) without analyzing the relation to a full constraint leaves a large field to 
be explored, which can give valuable consistency checks. We will not do this systematically 
in this paper but rather discuss a mathematical issue that arises in any such case: initially, 
one has to deal with difference equations of non-constant step-size which can be treated 
either directly or by tranforming a non-equidistant difference equation to an equidistant 
one. We first illustrate this for ordinary difference equations since partial ones, as they 
arise in anisotropic models, can often be reduced to this case. 

3.1 Ordinary difference equations of varying step-size 

Let us assume that we have an ordinary difference equation for a function tp^, which appears 
in the equation with /i-dependent increments ip^+sj\fi(fi)-'^ ■ To transform this to a fixed step- 
size, we introduce a new variable fl{fi) such that /2(/i+(5/A/i(/i)) = fl{fi) + 6fi'/J\fi{fi) + 0{6^) 
has a constant linear term in 6. (For the isotropic equation, A/i is the vertex number only in 
one direction. The total number of vertices in a 3-dimensional lattice is given by A/" = ATf .) 
This is obviously satisfied if we choose /i(/i) := A/i(z/)dz/. We then have 

^t^+S/XiM = = ^/i+5+J]-2^5Wi''-i)/Ar,' (26) 

1 \f' ~ 

= %^s + -S'j^,^P' + 0{6') (27) 

■^This behavior is introduced independently in 35J where "effective" equations, obtained by replacing 
connection components in the classical constraint by sines and cosines of such components according to 
how they occur in the quantized constraint, are analyzed for the Schwarzschild interior. The results are 
complementary to and compatible with our stability analysis of the corresponding difference equations 
below. We thank Kevin Vandersloot for discussions on this issue. 
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where A/'i*'' denotes the i-th derivative of A/i. Thus, up to terms of order at least 5^ the new 
equation will be of constant step-size for the function ipp_ := (The derivative ip' by 

/i may not be defined for any solution to the difference equation. We write it in this form 
since such terms will be discussed below in the context of a continuum or semiclassical 
limit where derivatives would exist.) 

It is easy to see that, for refining lattices, the additional terms containing derivatives of 
the wave function are of higher order in h and thus correspond to quantum corrections. For 
7Vi(/i) oc /i'^ as a positive power of /i, which is the expected case from lattice refinements 
related to the increase in volume, we have 



relating /i to an isotropic triad component p = 4:7ijipfi/3 as it occurs in isotropic loop 
quantum gravity Moreover, 

^, dtp dfidip 1 dip i 1 „^ 
d/i dp, dfi A/i(/i) dyU 2A/i(/i) 

in terms of a curvature operator c = 87Ti'jGh/3d/dp = 2id/dfi which exists in a continuum 
limit fTT]. Thus, 

-r-F^ip OC - Cip 



A/? \P. 

With q positive (or just larger than —1/2) for a refining lattice, there is a positive power 
of showing that additional terms arising in the transformation are quantum corrections. 

This has two important implications. First, it shows that the correct classical limit 
is obtained if lattices are indeed refined, rather than coarsened, since q is restricted for 
corrections to appear in positive powers of h. In anisotropic models, as we will see, the 
behavior is more complicated due to the presence of several independent variables. An 
analysis of the semiclassical limit can then put strong restrictions on the behavior of lat- 
tices. Secondly, we can implicitly define a factor ordering of the original constraint giving 
rise to the non-equidistant difference equation by declaring that all quantum correction 
terms arising in the transformation above should cancel out with factor ordering terms. 
We then obtain a strictly equidistant equation in the new variable fx. For example, a 
function A/i(/i) = gives jl oc such that the transformed difference equation 

will be equidistant in volume rather than the densitized triad component. For this special 
case, factor orderings giving rise to a precisely equidistant difference equation have been 
constructed explicitly in [251 El] • 



3.2 Number of vertices proportional to transversal area 

A simple difference equation results if the number of vertices is proportional to the transver- 
sal area in any direction]^ In the x-direction we have transversal surfaces given by symmetry 



* Since this refers to the area, it is the case which agrees with the motivation of [25l [34] . 
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orbits of area pc, using the line element ([3]), and thus oc r. Transversal surfaces for 
an angular direction are spanned by the x- and one angular direction whose area is pf,, 
giving J\f^ oc /i. Each minisuperspace direction has a step-size which is not constant but 
independent of the other dimension. Moreover, due to the simple form one can transform 
the equation to constant step-size by using independent variables and fx^ instead of r 
and /i. Illustrating the general procedure given before, a function ■^pT^2,^J.^ acquires constant 
shifts under the basic steps, 

i'iT+nS/T)2,{/i+mS/ii)2 = ^T'2+2n5+n262/T2,fi'2+2mS+m'2S2/ii'2 

up to terms which can be ignored for large r and fi. This is sufficient for a straightforward 
analysis in asymptotic regimes. Moreover, higher order terms in the above equation come 
with higher derivatives of the wave function in the form 

since q = 1 compared to the discussion in Sec. 13. 1[ Due to the extra factors of h (or even 
higher powers in further terms in the Taylor expansion) any additional term adding to the 
constant shift of can be attributed to quantum corrections in a semi classical limit. 

Accordingly, such terms can be avoided altogether by a judicious choice of the initial factor 
ordering of operators. 



3.3 Number of vertices proportional to extension 

Geometrically more intuitive, and as we will see below dynamically more stable, is the 
case in which the number of vertices in each direction is proportional to the extension of 
that direction measured with the triad itself. This gives M-a oc \/[r| and oc ^/ \/\t\, 
using the classical co-triad (jl]). (One need not worry about the inverse r since the effective 
treatment of lattice refinements pursued here is not valid close to a classical singularity 
where an already small lattice with a few vertices changes. Singularities in general can 
only be discussed by a direct analysis of the resulting difference operators. Since only 
a few recurrence steps are necessary to probe the scheme around a classical singularity, 
equidistant difference operators are not essential in this regime. They are more useful in 
semiclassical regimes where one aims to probe long evolution times as in the examples 
below. Similar remarks apply to the horizon at /i = which, although a classical region 
for large mass parameters, presents a boundary to the homogeneous model used for the 
Schwarzschild interior.) The behavior is thus more complicated than in the first case since 
the step size of any of the two independent variables depends on the other variable, too. 
First, it is easy to see, as before with quadratic variables, that the volume label uo = ^\/\t\ 
changes (approximately) equidistantly with each iteration step which is not equidistant for 
the basic variables and r. But it is impossible to find a second, independent quantity 
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which does so, too. In fact, such a quantity /(/x, r) would have to solve two partial 
differential equations in order to ensure that 



changes only by a constant independent of r and /i. This implies d^f{fi,T) oc a/|t| and 
r) oc /i/ a/pi whose only solution is /(/i, r) oc /Ua/PI which is the volume uj. 
We thus have to deal with non-equidistant partial difference equations in this case which 
in general can be complicated. A possible procedure to avoid this is to split the iteration 
in two steps since an ordinary difference equation can always be made equidistant as 
above (cancelling quantum corrections by re-ordering). We first transform r to the volume 
variable uj which gives, up to quantum corrections, constant iteration steps for this variable. 
With the second variable still present, a higher order difference equation 

-C_(/i,u;V/^^)V^Mi+25M,'^+<5 ~ 2Co(/i,t^7Ai^)(l + 27^5 V/^^^)/^^'^,^ 
+Co(/i,a;VAi^)(l - 25 /uj)nip^^i_^sM,Lu-45 = (28) 

results with 

Co(/i,u;V/^^) = + (29) 




C±(/x,^V/x2) = 26\^l + ^l±^j (30) 

derived from the original coefficients (12^ . The structure of this difference equation is quite 
different from the original one: not only is it of higher order, but now only one value of the 
wave function appears at each level of u, rather than combinations of values at different 
values of fi. Note also that only the coefficient of the unshifted ^/'^^i^ depends on /i. This 
form of the difference equation is, however, a consequence of the additional rotational 
symmetry and is not realized in this form for fully anisotropic Bianchi models as we will 
see below. 

Proceeding with this specific case, we have to look at wave functions evaluated at 
shifted positions /i(l + m6/uj) with integer m. At fixed u = uq, we are thus evaluating 
the wave function at values of fi multiplied with a constant, instead of being shifted by a 
constant as in an equidistant difference equation. This suggests to use the logarithm of 
fi instead of fi itself as an independent variable, which is indeed the result of the general 
procedure. After having transformed from t to u already, we have to use r as a function 
of fi and u in the vertex number Af^, which is T{fi,Lj) = {oj/fi)'^ after using u = fiy/r. 
Thus, J\f^{fi,T{fi,uj)) = a/ T{fi, uj) = uj/fi now is not a positive power of the independent 
variable fi and we will have to be more careful in the interpretation of correction terms 
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after performing the transformation. (The lattice is coarsened with increasing anisotropy 
at constant volume.) Naively applying the results of Sec. l3.1l to q = —1 would suggest that 
corrections come with inverse powers of h which would certainly be damaging for the correct 
classical limit. However, the factors change due to the presence of the additional variable 
uq even though it is treated as a constant. We have Mff/^f^ = -I/ujq = -(74/2)3/2/\/o 
in terms of the dimensionful volume V, while it would just be a constant —1 without the 
presence of cu. The additional factor of ^^/^ ensures that corrections come with positive 
powers of h for the correct classical limit to be realized. 

For any Uq, we thus transform 4'^(i+mS/ujo) equidistant form by using ip^ = ipf_i(fi) with 
= log/i. This transformation is possible since the second label Uq is now treated as a 
constant, rather than an independent variable of a partial difference equation. (Recall that 
for the type of difference equation discussed here there is only one variable, the volume, 
which is equidistant under all of the original discrete steps.) Despite of negative powers of 
some variables in the vertex numbers, we have the correct classical limit in the presence 
of a;. As before, the transformation is exact up to higher order terms which are quantum 
and higher order curvature corrections. Defining the original constraint operator ordering 
implicitly by the requirement that all those terms are cancelled allows us to work with an 
equidistant difference equation. 



3.4 Bianchi models 

As mentioned before, the transformed difference equation does not become higher order 
for fully anisotropic Bianchi models. In this case, we have three independent flux labels 
/i/, / = 1,2,3, and vertex numbers A/}. Using vertex numbers proportional to the spatial 
extensions for each direction gives A/i = a/ /i2/^3/A''i, A/2 = \/ Hi 1^3/ 1^2 and As = a/ Aii/i2//^3- 
As in the difference equation for the Schwarzschild interior, the difference equation for 
Bianchi models [22] uses values of the wave function of the form ^/'/^i+2<5/M,/i2+2<5/M,/i3- 
can again see easily that the volume lu = \/\jMj^2fJ^\ behaves equidistantly under the 
increments, 

cu(^i + 25/Ari,^2 + 2(5/Ar2,/i3) = ^ (^fii + 25^^) (^f^2 + 2'^^^) -"3 

= ^J^Jil^^2^^3 + 4:6^fiifi2fi3 + 4:6^ = uj + 26 + 0{6^) . 
The leading order term of the difference equation in uj results from a combination 

,/i2+2<5/A'2,M3+2<5/A/'3 + C2ll'fMi+25/Afi,fi2,t^3+25/Afs + C'^1pf^^^25/Afi,fi2+25/Af2,l^3 
i,tM2+2S/J\f2,i^+2S + C2^fj,i+25/Afi,fi2,uj+25 + C3'?/'^i+25/M,/i2+2<5/A^2,'^+25 
= Cl'^^ti,/i2(l+2<5/'^),w+25 + C2'^^ti(l+2<5/w),At2,a;+2<5 + C^lpf^ 

1 (l+2(5/w),^2 {1+26 / Lij) ,uj+25 

= : C+?/;a;+25(Atl,/i2) 

where we used 1/A/i = A//ii//i2yU3 = /^i/"^ and defined the operator (7+ acting on the 
dependence of ip on fii and fi2- Thus, unlike for the Schwarzschild interior the difference 
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equation does not become higher order in ui, and the highest order term does have a 
difference operator coefficient in the remaining independent variables. 

The recurrence proceeds as follows: We have a partial difference equation of the form 

with difference operators C± and Cq acting on the dependence on /^i and /i2- In terms of 
initial data at two slices of uj we can solve recursively for Coipujif^i, (^2) +C'-^a;-2<5(/^i, /^2) =: 
112) and then, in each w-step, use boundary conditions to solve the ordinary difference 
equation 

Although the operator (7+ itself is not equidistant, this remaining ordinary difference 
equation can be transformed to an equidistant one by transforming /ii and ^2 as in 
Sec. 13.11 (using that uj is constant and fixed for this equation at any recursion step). With 
fi3{fii, fi2,uj) = a;^//ii/i2, we have lattice spacings A/i(/ii, /i25 1^) = a;//ii and A/'2(/ii, /i25 1^) = 
Lj/lJ,2 in terms of 00 which are already independent of each other. The two remaining vari- 
ables fii and fi2 are thus transformed to equidistant ones by taking their logarithms as 
encountered before. 

Note the resemblance of the new variables, volume and two logarithms a metric compo- 
nents at constant volume, to Misner variables [36]. This observation may be of interest in 
comparisons with Wheeler-DeWitt quantizations where Misner variables have often been 
used, making the Wheeler-DeWitt equation hyperbolic. 

4 Application: Stability of the Schwarzschild interior 

Now that we have several possibilities for the lattice spacings, we consider their effect on 
the solutions of the Hamiltonian constraint. In particular, these solutions may have un- 
desirable properties reminiscent of numerical instabilities, as it was indeed noticed for the 
original quantization of the Schwarzschild interior in [|28J. Also problems in the presence 
of a positive cosmological constant, described in the introduction, are of this type. Re- 
call that when one wishes to solve an ordinary differential equation, for example, there 
are various discrete schemes that ensure errors do not propagate as the number of time 
steps increases. Here we are in the opposite situation - instead of having the freedom to 
pick the discrete version of a continuous equation, the discrete equation itself is what is 
fundamental. Thus, like a badly chosen numerical recipe, some choices of the functions A^t 
and in the constraint equation may quickly lead to solutions that are out of control, 
and increase without bound. To test for this, we will use a von Neumann stability analy- 
sis [28j on the possible recursion relations. The essential idea is to treat one of the relation 
parameters as an evolution parameter, and decompose the rest in terms of orthogonal 
functions, representing "spatial" modes of the solution. This will give rise to a matrix that 
defines the evolution of the solution; if the matrix eigenvalues are greater than unity for 
a particular mode, that mode is unstable. In particular, a relation J2'k=-M '^n+ki^n+k = 
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is equivalent to a vector equation of the form Vn = Q{n)vn-i, where the column vector 
Vn = ii^n+M, i^n+M-i, ' ' ' 5 4'n-M+i)'^ ■ The cvolutiou of an eigenvector w of the matrix Q{n) 
is given by w„ = XwWn-i- Thus, when the size of the corresponding eigenvalue |A^| > 1, 
the values in the sequence associated to w will grow as well. 

With this in mind, we consider the choices of and A^^ discussed previously, starting 
with the case = t and = fi. In the large n, r limit for this choice, the coefficients of 
the Hamiltonian constraint become 



C±(/x,r) 



T ~ 



-3/2 ' 



In the asymptotic limit, the coefficients of the 4'^±25/^i,T and ip^^r terms go to Co(/U, r)/i. 
As we saw in Section 3.2, we can choose a different set of variables in which the step sizes 
are constant (up to ordering of the operators). Plugging these asymptotic values into the 
Hamiltonian constraint, and changing variables to fl = and f = gives 

4:f{lljfi+25,f+25 - ^il-25,f+2S + V^/i-2<5,f-25 - ^fL+2S,f-2s) + fl{lpjl+45,f - '2lpjl,f + ^A-4<5,^) = 0. 

Because all the step sizes now are constants depending on 6, we define new parameters 
m, n such that ft = 2m5 and f = 2n5. Using m as our evolution parameter and n as the 
"spatial" direction, we decompose the sequence as 'il^2m.5,2n5 = Umexp(mu;). With this new 
function, the recursion relation is written as 



2in{Un+l - Un-l] 

This is equivalent to the vector equation 



m sm f 



0. 



" Un+l ' 




- Un . 





-f^sinO 1 

2n 

1 







= Q{m,n) 


Un 




- Un-l . 


■ Un-l - 



(31) 



The eigenvalues of the matrix Q are 



Ah 



-im sin 6 ± y IQ-n? — iv? sin^ 9 
An 



When the discriminant IQn^ — w? siv? 6 > 0, then |A| = 1, and the solution is stable; 
however, there are unstable modes when 16n^ — m^sin^^ < 0. The most unstable mode 
corresponds to the choice sin^ = 1, giving instabilities in terms of the original variables 
when /i > 2r. In this regime, all solutions behave exponentially rather than oscillating. 
This region includes parts of the classical solutions for the Schwarzschild interior even for 
values of /i and r for which one expects classical behavior to be valid. The presence of 
instabilities implies, irrespective of the physical inner product, that quantum solutions 
in those regions cannot be wave packets following the classical trajectory, and the correct 
classical limit is not guaranteed for this quantization, which is analogous to that introduced 
in [251131]. 
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The situation is different when we consider the choices = y |t| and A^^ = jj,/ \/\t\^ 
where we will find a lack of instability. There is no choice of variables that allows us to 
asymptotically approach a constant spacing recursion relation, because of the mixing of 
the /i and r variables in the step size functions. Thus, we will make the assumption that 
in the large /x, r limit, the solution does not change much under step sizes 5N~^ and 5N^^. 
To see how this affects the resulting stability of the solutions, we will look at a simpler 
example first. If we start with the Fibonacci relation R^- = tpr+i ~ '4't ~ "^t-i = 0, then 
the two independent solutions are of the form ip^ = , where k is either the golden ratio 
= (1 + \fE)/2 or else — Only the latter solution meets the criterion for stability, 
since |0| > 1. When we change this relation to 

Rr = V'r+l/r" " V^r " V'r-l/r" = 0, (32) 

with n 7^ 1, the situation changes - only one of the two solutions outlined above will solve 
the relation asymptotically. In particular, when we examine the error R^ we get when we 
plug kT' into the altered relation (!32|) . i.e. 

the error is proportional to ipr itself. As r — > oo, therefore, the error for the k = (p solution 
grows without bound, while that of k, = — goes to zero. Thus, we see in this situation 
a relation between the stability and the asymptotic behavior of a solution. 

Returning to the Schwarzschild relation, in the large /x, r limit the coefficient functions 
of the recursion relation are to leading order 

C±(/i,r)~45, Co(/i,r)~-. 



In turn, the relation itself becomes 



+ (^/.+45/v/7,r - 2V^^,r + V^^-45/v/7,t) = 0. 

From this point on, we assume that we have a solution to this relation which does not vary 
greatly when, for example, /i is changed by ±26/y/JI, and similarly for r. Both A^^. and A^^ 
are constant to first order in shifts fi ± 25 and similarly for r, in the asymptotic limit. 
Thus, we assume that a = 26 and P = 26 N^^ are constants, and use the scalings 
/i = am and r = f3n. When this is done, we get an equation similar to the case when 
Nx = T and A^^ = fi, but with constant coefficients; this is the crucial difference that allows 
stable solutions to the case here. Using the decomposition ipam,/3n = Un exp{im9), we arrive 
at the matrix equation 



' Un+l ' 




- Un . 





I sin 6 1 1 r Ur. 



1 J 



Un-l 



(33) 



The matrix here has eigenvalues A with |A| = 1 for all m, ra, so the solution is stable. Using 
arguments as in the Fibonacci example, the non-equidistant equation of the second scheme 
is shown to be stable. 



17 



5 Conclusions 



Following [9j, we explicitly introduced loop quantum cosmological models which take into 
account the full lattice structure of inhomogeneous states. Such lattices are in general 
refined by adding new vertices when acting with the Hamiltonian constraint. Thus, also 
dynamical equations even in homogeneous models should respect this property. Several 
interesting features arose: One obtains non-equidistant difference equations which, when 
imposed for functions on the whole real line as in isotropic loop quantum cosmology, 
are more restrictive than equidistant ones due to the absence of superselected sectors. 
This leaves the singularity issue unchanged since for this one only needs to consider a 
few steps in the equation. But a stability analysis of solutions and the verification of the 
correct classical limit in all semiclassical regimes can be more challenging. We presented an 
example for such an analysis, but also introduced a procedure by which one can transform 
the resulting equations to equidistant ones up to quantum corrections, which is sufficient 
for a semiclassical analysis. Interestingly, properties of the transformation itself provide 
hints to the correct semiclassical behavior. As a side-result, we demonstrated that one 
particular version of lattice refinements naturally gives rise to Misner-type variables. 

It is our understanding that this general procedure of defining lattice refining models 
mostly agrees with the intuition used specifically in isotropic models in [2B], and adapted 
to anisotropic ones in [51] Jf] However, there are some departures from what is assumed in 
[25] . First, we do not see indications to refer to the area operator while the area spectrum 
was not only used in [25| to fix the constant 6 and the volume dependence of the step size 
but in fact provided the main motivation. Secondly, due to this motivation [25] presents 
a more narrow focus which from our viewpoint corresponds to only one single refinement 
model. It has a vertex number proportional to volume, which is a limiting case not realized 
by known full Hamiltonian constraints, and puts special emphasis on geometrical areas to 
determine the vertex number. Finally, commutators for inverse volume operators are to be 
treated differently from [25], taking into account a lattice refining model which would not 
be possible in a purely homogeneous formulation. As shown in the appendix, this enlarges 
expected quantum corrections to the classical functions. 

We have discussed similar cases for illustration here, but keep a more general view- 
point on the refinement as a function of volume. A preliminary stability analysis for 
the Schwarzschild interior, consistent with [35] indeed suggests that a behavior different 
from what is suggested in [25j is preferred, which indicates that models can provide tight 
conditions for the general analysis of quantum dynamics. We emphasize that stability 
arguments as used here are independent of physical inner product issues since they refer 
to properties of general solutions. A general analysis as started here allows detailed tests 
of the full dynamics in manageable settings, which can verify the self-consistency of the 
framework of loop quantum gravity — or possibly point to limitations which need to be 
better understood. 

^We thank A. Ashtekar for discussions of this point. 
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A Inverse volume terms in homogeneous models and 
lattice refinement 

We have seen in this paper, following [9] , that Hamiltonian constraint operators with triad- 
dependent parameters in holonomies allow one to model lattice refinements faithfully, with 
interesting results and some improvements over the original non-refining models. However, 
as always there are also some features of inhomogeneous states and operators which are 
not present and difficult to mimic in homogeneous models. Thus, even models generalized 
in this way by allowing for lattice refinement effects have to be interpreted with great care. 
While qualitative effects can be investigated fruitfully to test the full framework, there is 
no basis for drawing quantitative conclusions. The prime example is that of commutator 
terms which appear ubiquitously in composite operators of loop quantum gravity, such as 
the coefficients of difference equations or also matter Hamiltonians. 

In the main construction of this paper we used holonomies associated with links of a lat- 
tice, rather than edges of a fixed coordinate length. This allows us, effectively, to take into 
account lattice refinements which change the number of vertices. It applies to holonomies 
( |T6|) along a closed loop used to quantize curvature components which determine the step 
size of difference equations, but also to the link holonomies used in commutators to quan- 
tize inverse triad components based on (fT^ . What is not modeled in homogeneous models 
is the fact that a lattice operator makes use of the local volume K, at a given vertex v 
where the commutator is acting, rather than the total volume V = Ylv of the whole 
box in which the lattice is embedded. In a fully inhomogeneous setting the difference 
does not matter since volume contributions from vertices not touched by the edge used in 
a commutator drop out in the end, [/^e, ^] = J^veel^e^^v]- But in homogeneous models 
there is a difference since volume contributions from different vertices, in an exactly ho- 
mogeneous setting, are all identical. Thus, the total volume V = MV^ is the number of 
vertices multiplied with the local volume V^. Then, [h, V^] rather than [h, V] is expected as 
the contribution to constraint operators from the inhomogeneous perspective. In homoge- 
neous models as in [25|, on the other hand, [/i, V] is more straightforward to use. We now 
show that without corrections this would imply crucial deviations from the inhomogeneous 
behavior. 

It is easy to see that commutators differ depending on whether the local or total volume 
is used. For simplicity of the argument, we proceed with an isotropic situation where 
V = Ipp/^ in terms of the basic isotropic densitized triad component p. A local lattice flux, 
for a surface 5* intersecting only a single link, would be p = J^(Pyp = i^p/Ll = p/N"^^^ 
for links of coordinate length such that = = \p\'^^'^/J\f is the local volume. 
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(We again use a coordinate box of size Lq and introduce rescaled flux variables p = L^p.) 
Isotropic states are spanned by e'^'^/^ where /i G M is related to the triad eigenvalues by 
= 47r7£p/i/3. Using a link holonomy h ~ e''^oc/2 _ gic/2A/'i/3 ^j-^j^j-^ multiplication 
operator increases ^ by l/M^^^ , a commutator with the total volume will have eigenvalues 
of the form 

h-\h,V] ~ \/(/i + l/Ari/3)_y(^_i/Ari/3) = |p+47r74/3Ari/=^|3/2_|^_4^y2/3_^i/3|3/2_ 

(34) 

If the local volume is used, on the other hand, we have to refer to local edge labels /Xg = 
3p/47r7£p, rather than using the total p. Thus, 

h'Ahe.V,] ~ K(/^e + l)-K(/^e-l) = |p + 47r74/3|3/2-|p-47r74/3|3/2 

= M-\\p + 47r74Ar2/V3|='/2 _ |p _ 4;r74Ar2/V3|3/2) (35) 

For large volume, p ^ A/", both expressions give the correct classical limit ^/\p\ 
expected from |e«c/2Ari/3^ |p|3/2j^ However, quantum corrections, i.e. deviations from this 
classical limit for finite p, are much larger for the second version using the local volume as 
it would occur in an inhomogeneous quantization. The smooth classical function dV/ dp in 
a Poisson bracket appears in discretized form by the large step-size N"^/"^ in (135!) rather than 
the small one N~^^^ in flMl) . Perturbative corrections, derived by Taylor expanding the 
difference terms and keeping higher order corrections to the classical expression, are thus 
larger. (This can have cosmo logical implications [371 EH [39] •) Non-perturbative effects as 
observed for the inverse scale factor operator in isotropic loop quantum cosmology which 
has an upper bound at finite volume |10], start to arise for p ~ M"^^^ when the local 
volume is used but only at the much smaller p ~ Af~^^^ for the total volume. Since only 
the local volume is relevant for inhomogeneous quantizations, quantum corrections from 
inverse volume operators can be large. 

Unfortunately, this effect is more difficult to mimic in exact homogeneous models unlike 
the behavior of holonomies under lattice refinements and has therefore been overlooked in 
|25] . The connection components appearing in holonomies can simply be divided by a 
function A/" of triad components to implement shrinking edges due to subdivision. This is 
not possible for the volume itself to use a local version in a homogeneous model since, if 
we would divide the total volume by the appropriate function of triad components, only 
a constant would remain for an Af proportional to volume and the commutator would be 
zero. The only way to have this effect faithfully implemented in a homogeneous model is 
to use higher SU(2) representations for holonomies in commutators but not for holonomies 
used in the loop to quantize curvature components. (This is not possible if one writes the 
constraint as a single trace, tT{hahe[h~^ ,V]) but can easily be done using the equivalent 
form tT{Tiha)tT{Tihe[h~^,V]). We emphasize that higher representations for commutators 
are advocated here only in exactly homogeneous models to mimic inhomogeneous effects. 
Fully inhomogeneous operators usually need not refer to higher representations.) In a rep- 
resentation of spin j, matrix elements of holonomies contain exponentials exp(imc/2A/'^/^) 
with —j<m<j, which increases the shifts in volume labels resulting from commutators. 
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Resulting expressions for commutators can be found in [HlllS]. If the representation label 
j is of the order J\f, effects as they result from lattice refinements and using the local vol- 
ume are correctly implemented. Accordingly, corrections from inverse triad components 
quantized through commutators are much larger than they would otherwise bej§ 
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